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I consider the effects of enforced dimerization on random Heisenberg antiferromagnetic S = 1 
chains. I argue for the existence of novel Griffiths phases characterized by two independent dynamical 
exponents that vary continuously in these phases; one of the exponents controls the density of 
spin-1/2 degrees of freedom in the low-energy effective Hamiltonian, while the other controls the 
corresponding density of spin-1 degrees of freedom. Moreover, in one of these Griffiths phases, the 
system has very different low temperature behavior in two different parts of the phase which are 
separated from each other by a sharply defined crossover line; on one side of this crossover line, the 
system 'looks' like a S = 1 chain at low energies, while on the other side, it is best thought of as 
a S — 1/2 chain. A strong-disorder RG analysis makes it possible to analytically obtain detailed 
information about the low temperature behavior of physical observables such as the susceptibility 
and the specific heat, as well as identify an experimentally accessible signature of this novel crossover. 



I. INTRODUCTION 

Quantum spin chain systems are known to exhibit 
many interesting states of matter at low temperature 
that arise from the fascinating interplay between quan- 
tum fluctuations, correlation effects, and the effects of 
quenched disorder. They are particularly interesting 
from a theoretical point of view, as they provide ex- 
perimentally realizable examples in which this interplay, 
common to many other condensed-matter systems, can 
be studied in detail. 

Among the most studied such systems is the Heisen- 
berg antiferromagnetic (HAF) S = 1 chain (for an exper- 
imental review see Ref. [l]) with Hamiltonian: 

H = J^Si-S i+1 , (1) 

i 

where J is positive, and the operators S represent the 
spin angular momentum of spin-1 objects on site i of a 
one-dimensional lattice. As Haldanecl showed many years 
ago, quantum fluctuations in the S = 1 HAF chain are 
strong enough to rule out even quasi-long range antifer- 
romagnetic order of the type present in S = 1/2 chains, 
and the system is in the so-called 'Haldane' phase charac- 
terized by a singlet ground-state and gap to all bulk ex- 
citations. However, this phase possesses a kind of 'topo- 
logical' ordem that results in the presence of almost free 
sub-gap boundary spin-1/2 degrees of freedom in a large 
but finite system of length L with free ends; these spin- 
l/2s are coupled to each other by an effective exchange 
coupling J c ff ~ (— l) L+1 e~ cL (with L in lattice units, 
and c proportional to the inverse ground-state correla- 
tion length), and will play a crucial role in much of the 
analysis in this paper. 

A very intuitive and useful picture of the Haldane 
phase is the 'valence-bond' solida description: The ba- 
sic idea is to think of each spin-1 as a (symmetrized) pair 
of spin-1/2 objects; each site thus has two 'indistinguish- 
able' spin-1/2 objects on it. A good description of the 



ground state is then obtained by pairing one of the spin- 
1/2 at any site into a singlet state with a spin-1/2 on its 
neighbor to the right, and pairing the other spin-1/2 into 
a singlet state with a spin-1/2 on its neighbor to the left. 
This gives a state with one singlet bond across each link 
of the lattice, i.e the (1,1) state. The topological order 
present in the Haldane phase can now be thought of sim- 
ply as the ability to 'walk' across the length of the system 
along singlet bonds without encountering any breaks, and 
the unpaired spin-1 /2s at each end of a chain with free 
boundaries model the sub-gap end-states characteristic 
of the Haldane phase. 

The stability of this phase to various experimentally 
relevant perturbations has also been investigated. For 
instance, it is knownH that the Haldane phase is stable to 
weak enforced dimerization 5 in the exchange couplings 
(5 parametrizes the extent to which the even bonds are 
stronger than the odd bonds), and there are no qualita- 
tive differences in the low-energy properties of the system 
until \S\ exceeds a critical value S c . Beyond this point, the 
ground state changes character, and the system enters a 
different gapped state without the topological order or 
the subgap spin-l/2s on free ends. For positive 5 > S c , 
this dimerized phase corresponds to the (2,0) valence- 
bond state in which both spin-1/2 degrees of freedom at 
any even-numbered physical site form singlets with spin- 
l/2s on its neighbor to its right, and the number of singlet 
bonds across each link of the lattice alternates between 
two and zero (similarly, the (0,2) state is favored by the 
system when 8 < —5 C ). 

The effects of quenched randomness R in the ex- 
change .couplings (with all J, still positive) have also been 
studiedoQ at 8 — (in the presence of randomness, 5 is 
defined in terms of the average values for even and odd 
bonds). The topological order characteristic of the Hal- 
dane state persists for R less than a critical value R c , and 
the Haldane phase is thus stable to weak randomness. On 
the other hand, a different, randomness dominated state 
is stabilized for R > R c . In this 'Spin-1 Random Singlet' 
(RSi) state, the interplay of disorder and quantum me- 
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FIG. 2. A cartoon for Griffiths effects in the Gapless Hal- 
dane phase at 8 — 
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FIG. 1. Schematic phase diagram as a function of ran- 
domess (R) and dimerization (5). Details in text 



chanics locks each spin-1 into a singlet state with some 
other spin-1; the two spins in a given singlet pair can 
have arbitrarily large spatial separation, with disorder 
determining the actual choice of partner for each spin;u 
in the valence bond picture, this large R state is thus 
characterized by a particular random pattern of double 
bonds. 

Although the topological order survives for all R < R c , 
it turns out that disorder has a dramatic effect on the 
low-energy properties near the transition, resulting in ar- 
bitrarily low-energy excitations with a power-law density 
of states controlled by a single non-universal (R depen- 
dent) dynamical exponent z. These low-energy excita- 
tions are associated with large anomalous regions of the 
sample in which the couplings locally favor a pattern of 
singlet bonds more characteristic of the other 'nearby' 
phase(s) (see below). Such 'Griffiths effects', whereby 
rare disorder-induced fluctuations in the interactions over 
extended regions of space give rise to non-universal sin- 
gular contributions that the dominate low-energy prop- 
erties, are among the more interesting and ubiquitous 
aspects of the physics of low-dimensional random quan- 
tum systems, and below, I focus on precisely this physics 
for the general S ^ case. 



II. PHYSICAL PICTURE AND MOTIVATION 

Let us start with the overall structure of the phase di- 
agram in the (R,5) plane: To begin with, note that the 
nature of the ground state will not change qualitatively 
from that at S = as long as R < R c and 5 is small 
enough. On the other hand, by analogy with the parre- 
sponding situation in r andom HAF 5 = 1/2 chainsjiil the 
RSi state for R > R c will be unstable to adding a tiny 
amount of dimerization. For small 6 > 0, the resulting 



pattern of singlet bonds in the ground state will be the 
random analog of the (2,0) state of the pure system, and 
topologically indistinguishable from the latter. Likewise, 
the state for infinitesimal <5 < will be topologically iden- 
tical to the (0,2) state. Putting all this together, we are 
led to a phase diagram that looks something like the one 
sketched in Fig 0. The critical pointfitj at R = R c , 5 = 
is thus really a multicritical point at which the (2,0), (0,2) 
and (1,1) phases meet, while the RSi line is seen to be 
the critical phase boundary between the (0,2) and (2,0) 
phases. 

To complete this part of our discussion, we need to 
identify the universality class of the phase-boundary be- 
tween the (1,1) phase and either of the dimerized phases. 
To this end, note that the physics asymptotically close 
to ±<5 C at R — is that of a spin-1/ 2 chain with a small 
value of enforced dimerization: heuristically, this may be 
understood by first putting down a single valence bond 
on each even link (for 5 > 0) to partially 'screen' out the 
imbalance in the odd and even exchange couplings; this 
leaves behind .a 5=1/2 HAF chain with 5 e s ~ 6 — 6 C . 
Now, we knowa that adding disorder to the 5=1/2 HAF 
at S e s = leads to a spin- 1/2 Random Singlet (RSi/ 2 ) 
state (with a random pattern of single valence bonds sta- 
tistically identical to the pattern of double bonds in the 
RSi phase), and thus, the (5^0 phase boundaries be- 
tween the (1,1) state and the (2,0) or (0,2) states are lines 
along which the system is in the RS1/2 critical state. 

Of course, the topological labels of the different phases 
and the universality classes of various transitions are not 
the whole story, and we need to consider the role of Grif- 
fiths effects in various regimes to really understand the 
low-energy behavior of the phases. Consider, for starters, 
the Gapless Haldane regime at S = 0. The gapless spec- 
trum obtained in Refs |]j7| can be understood by thinking 
about the effect of a single rare region of length L in 
which the exchange couplings are such that the preferred 
pattern of singlet bonds in this region is more character- 
istic of a system at a nearby point in either the random 
(2,0) or (0,2) phases, or on the critical RSi line separating 
them (outside of this rare segment of length L, the sys- 
tem is in the 'typical' (IT) state). In all such cases, the 
central segment (see Fig ||) can be thought of as an 'insu- 
lating' barrier that separates two essentially semi-infinite 
chains in the Haldane phase, implying the presence of 
low-energy spin-1/2 degrees of freedom localized at the 
ends of the (1,1) segments. 

These spin-1 /2s are coupled to each other across the 
'barrier' by an antiferromagnetic exchange coupling of 
order e~ ClL , while the probability for such an anoma- 
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FIG. 3. A caricature of Griffiths effects in the (1,1) phase 
at positive S 
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FIG. 4. Cartoon for Griffiths effects in the (2,0) phase 



lous disorder configuration to occur is also exponen- 
tially small: pl ~ e~ C2i . Averaging over different 
possibilities, with anomalous regions of varying length 
L > — ln(0)/ci, then yields a power-law for the density 
of spin- 1/2 degrees of freedom that model the physics be- 
low the energy scale Q: n±/2 ~ fi 1 / 2 , with a non-universal 
exponent z (here, and henceforth, f2 denotes some suit- 
able low energy scale much smaller than the microscopic 
scale Jtyp). Now, if two such anomalous regions occur es- 
sentially 'next' to each other [i.e separated by a segment 
with typical exchange couplings of length L < — ln(f2)), 
and if L is even, then the end spin-1/2 degrees of freedom 
of this intervening typical segment will bind ferromagnet- 
ically to each other producing an effective spin-1 as far as 
the low-energy dynamics below the cutoff scale f2 is con- 
cerned. The probability for this to happen is dominated 
by the requirement that two anomalous regions, that are 
usually separated by lengths of order SI -1 / 2 , occur essen- 
tially next to each other — this immediately implies that 
the density of spin-1 degrees of freedom in the effective 
Hamiltonian will scale as the square of the density of 
spin-1/2 objects, consistent (apart from a log-correction) 
with the RG results of Ref ||. 

It turns out that similar heuristic arguments can be 
employed to qualitatively understand Griffiths effects at 
S > in the (1,1) phase, i.e, closer to the phase boundary 
to the (2,0) state (see fig g): In this case, the dominant 
disorder-induced fluctuation will consist of large regions 
of the sample that locally want to be in the (2,0) state — 
the probability of their occurrence will control the ex- 
ponent z 1 / 2 that determines the density n 1 / 2 ~ SI 1 / 21 / 2 
of low-energy spin-1/2 degrees of freedom below scale 
f2, while much rarer regions that locally prefer the (0,2) 
phase or the RSi critical state will only provide sub- 
dominant corrections as far as the value of Zi/% is con- 
cerned. 

However, it is clearly impossible (simply for parity rea- 
sons) to have two anomalous segments in the (2,0) phase 
separated from each other by an even segment — one of 
these anomalous regions necessarily has to be of a much 
rarer sub-dominant type (either a (0,2) region, or one 
with a spin-1 random singlet pattern of valence bonds) 
in such a situation. These are precisely the configurations 
that behave at low energies as a spin-1 object, and thus, 
the spatial density of spin-1 degrees of freedom below 



scale fl vanishes faster than Q 2 / Zl / 2 . Introducing a sec- 
ond exponent z\ < 21/2/2 to describe this behavior, we 
therefore arrive at n\ ~ J7 1 / 21 . Moreover, it is also clear 
from this argument that z\ 1% will diverge as we approach 
the transition at non-zero 8 into either the (2,0) or (0,2) 
phases (after all, what was the dominant rare fluctua- 
tion becomes the typical configuration as one crosses the 
transition!), while z\, whose value near the RS1/2 phase 
boundary is primarily determined by the sub- dominant 
fluctuation, will remain finite at this transition. 

In this analysis of the (1,1) phase, I have been un- 
avoidably led to a description in terms of two independent 
dynamical exponents. This is, of course, formally quite 
interesting — however, since Z\ < £1/2/2 throughout this 
phase, this additional 'tuning-knob' does not lead to any 
qualitative changes in the low-energy physics as we move 
around in this phase. On the other had, the situation 
in the (2,0) (or (0,2)) phase in the vicinity of its phase 
boundaries is dramatically different, and this is what I 
turn to next. 

To begin with, let us focus on two representative points 
in this phase (see Fig ^) , lying on an arc drawn around 
the multicritical point and going from the RSi critical 
line to the RS1/2 critical line. Consider first the point Pi 
very close to the RSi line: The low-energy dynamics at 
this point will be dominated by rare large regions that 
locally want to be at a nearby point in the (0,2) phase, 
embedded in a more typical background which is in the 
(2,0) phase (see fig||). Such an anomalous region clearly 
has residual spin-1 objects at either end, and the low- 
energy dynamics below scale fl will be dominated by a 
density n\ ~ fi 1 / 21 of these, with Z\ being controlled 
by the probability for such anomalies to occur at the 
corresponding length scales. 

Of course, there will also be sub-dominant disorder 
fluctuations that involve much rarer regions locally in 
the (1,1) phase; these have spin-1/2 degrees of freedom 
at their ends, and the low-energy dynamics will there- 
fore also have sub-dominant contributions from a density 
ni/2 ~ SI 1 / 21 / 2 of these (with Z\i% *C Z\ being controlled 
by the probability of occurrence for the (1,1) anomalies). 
Simply put, in this regime of the (2,0) phase, the low- 
energy behavior is more or less that of a spin-1 system. 
Next, consider the point P 2 chosen very close to the RS1/2 
phase boundary. The situation is now completely re- 
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versed, and the dominant Griffiths effects arise from rare 
regions locally at a nearby point in the (1,1) phase. Thus, 
Z1/2 will be much larger than z\ in this regime, and the 
low-energy behavior is predominantly that of a system 
made up of spin- 1/2 objects. 

Clearly, as we move along the arc from P\ to P2, the 
low-energy response will smoothly interpolate between 
these two extremes. Thus, the system in either the (2,0) 
or (0,2) phase has a curious ability to look qualitatively 
different at low energies in different parts of the phase, 
undergoing a metamorphosis from a spin-1/2 system into 
a spin-1 system as we move around in the phase. In 
particular, one expects that there will be some interme- 
diate point along the arc at which z\ = £1/2; at this 
point, both spin-1/2 and spin-1 degrees of freedom will 
play a role in determining the low-energy behavior of 
the system. [More generally, one expects a crossover 
line emanating from the multicritical point, along which 
m/2 ~ n% ~ SI 1 / 2 as Q — > 0.] 

These heuristic considerations clearly demonstrate 
that there is enough new and interesting physics at 5 
to warrant a more serious analysis, and this is what I 
turn to next. 



III. EFFECTIVE MODEL AND THE RG 
APPROACH 

To go beyond the qualitative arguments of the previ- 
ous section, one needs to analyze an appropriate effective 
model in a controlled manner. To fix the form of this ef- 
fective model, it is useful to specialize to a situation in 
which the bonds Ji in our Hamiltonian 



H — JiSi ■ 5,;+! 



(2) 



only take values 1, e s > 0, or e w > (s and w stand for 
'strong' and 'weak', with 1 > e s > e w ). For i even, Ji 
is 1 with probability 1 — p s (the Ji for i even are thus 
'strong'), and e s with probability p s , while, for i odd, Ji 
is 1 with probability 1 — p w (the Ji for odd i are thus 
'weak'), and e w with probability p w (with p w > p s ). For 
concreteness, it is also useful to take the p s / w and e s / w 
all much smaller than 1. In this case, it is appropriate to 
think in terms of a collection of fairly large segments of a 
pure S = 1 HAF chain in the Haldane state, coupled to 
each other by the e-bonds, and described in terms of pairs 
of the sub-gap boundary spin-1/2 degrees of freedom at 
their ends. 

Given that p w > p s , a typical segment of this kind 
will be, more likely than not, flanked on both sides by 
e w -bonds. Consequently, such a segment will be odd in 
length, and its low-energy behavior will be described by a 
pair of spin-l/2s coupled to each other by an 'odd' bond 
which is antiferromagnetic, and drawn from a distribu- 
tion that is calculable in terms of the length-distribution 
of the pure segments (the designation 'odd' has nothing 
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FIG. 5. Some examples of the different types of clusters 
that make up our effective model. 



to do with the parity of the segment length, and is chosen 
to conform to the notation of Ref. || when S — 0) — the 
bonds connecting this 'odd' pair to the rest of the system 
on either side are 'even' bonds of the w-type (the desig- 
nation 'even' is again chosen to conform with Ref [| when 
6 = 0; clearly all 'even' bonds are always antiferromag- 
netic). On the other hand, a segment flanked on one side 
by an e s -bond, and on the other by a e w bond, will be 
even in length, and its low-energy description will consist 
of a pair of spin-l/2s coupled by a ferromagnetic 'odd' 
bond (whose modulus has the same distribution as in the 
earlier antiferromagnetic case). Of course, this odd fer- 
romagnetic pair is flanked by one s-type even bond and 
another w-type even bond; odd ferromagnetic bonds thus 
have a s -flank and a w -flank. 

Moreover, and this is crucial, the other e-bonds in the 
vicinity of any such segment (that contributes a ferro- 
magnetically coupled odd pair of spin-l/2s to the effec- 
tive model) are again more likely to be of the w-type, and 
thus, there is an enhanced probability for finding another 
even-length pure segment (contributing another pair of 
spin-1/2 coupled by a ferromagnetic odd bond) close-by 
on the s-flank of this segment. In other words, a ferro- 
magnetic odd bond in the effective theory likes to have 
another ferromagnetic odd bond (or a spin-1 object; see 
below) close-by on its s-flank. Likewise, an odd-length 
segment flanked on both sides by e s bonds (described at 
low-energies by an odd antiferromagnetic pair of spin- 
l/2s flanked on both sides by s-type even bonds) is more 
likely than not to have nearby even-length segments in 
the Haldane phase; thus, in the effective theory, an odd 
antiferromagnetic bond flanked by s-type even bonds on 
both sides side likes to have odd ferromagnetic bonds (or 
spin-Is; see below) close by on both sides. Finally, two e 
bonds right next to each other will give rise to a spin-1 
object connected to the rest of the system by a ui-type 
even bond on one side, and a s-type even bond on the 
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other (obviously, such a spin-1 object also has a tendency 
to have another ferromagnetic odd bond (or spin-1) close- 
by on its s-flank). Clearly, one could equally well repre- 
sent this spin-1 by two spin-l/2s connected by a ferro- 
magnetic odd bond with strength much bigger than the 
cutoffJS'EI and in this sense, this is just a special case of 
a ferromagnetically coupled odd-pair — however, I prefer 
to explicitly introduce spin-1 objects into my description 
here. 

This dimerization induced clustering tendency of the 
ferromagnetically coupled odd pairs (or spin-Is) in the 
effective model immediately implies that individual pairs 
of spin-l/2s coupled by odd bonds (and their special case: 
spin-Is) are not the elementary constituents of our effec- 
tive theory (this should be contrasted with the S — case, 
where the low-energy effective model does consist of sta- 
tistical]* independent odd pairs coupled together by even 
bonds) Jj Clearly, the solution is to think not in terms of 
individual odd pairs, but in terms of clusters: The sim- 
plest cluster, the a-type cluster, is just a pair of spin-l/2s 
coupled by an antiferromagnetic odd bond, and flanked 
on both sides by w-type even bonds. Any other cluster of 
mass [i is made up of a string of [i antiferromagnetically 
coupled odd pairs of spin-1 /2s (// = , 1 , 2 , 3 . . .), each 
flanked by s-type even bonds on both sides (the total 
number of s-type even bonds is thus fj, + 1), terminat- 
ing at each end in a ferromagnetically coupled odd pair 
of spin-l/2s. Of course, either, or both, of the termi- 
nating pairs could equally well be replaced by a single 
spin-1 each, and there are thus three kinds of clusters: 
11 (with two spin-Is at the two ends), If (with a spin- 
1 at one end and a pair of spin-l/2s coupled with an 
odd ferromagnetic bond at the other end), and // (with 
ferromagnetically coupled odd pairs of spin-l/2s at both 
ends). |— | 

Our effective modeltil is thus made up of these four 
types of statistically independent clusters, connected to 
each other by intervening w-type even bonds. Further- 
more, the absolute values of all odd bonds, regardless of 
their sign, are expected to be governed by a single dis- 
tribution P . On the other hand, the w-type and s-type 
even bonds will have two different distributions, which 
we denote by P w and P s respectively. Furthermore, since 
each end of any cluster is expected to independently be 
either a spin-1, or a ferromagnetically coupled odd pair 
of spin-l/2s, the relative abundances of the 11, If, and 
// clusters can be parametrized by a single probability g 
for an end to be a spin-1 (with the probability for an end 
to be a ferromagnetically coupled pair of spin-l/2s being 
1 — g). Similarly, the density of the 'elementary' a-type 
clusters relative to the density of all other clusters can be 
parametrized by a probability f a for any given cluster to 
be of type a. Finally, since we are dealing with uncorre- 
cted disorder, it is reasonable to expect that the masses 
\x of the 11, If, and // clusters are all characterized by 
a single exponential distribution, with probability cx 
for any cluster to have mass /i. 

All these expectations are borne out by my detailed 



calculations below; for now, I only note that the effective 
model is completely specified by the three probabilities g, 
f a and q, the three probability distributions P s , P w and 
P a , and the number of w-type even bonds N w . Naturally, 
the initial values of these parameters and distributions 
(at the intermediate energy scale fio below which this 
model becomes applicable) are determined by the com- 
plicated interplay of dimerization and randomness in the 
microscopic model — fortunately, their precise values and 
functional forms will, for the most part, be unimportant 
as far as the physics at scales f2 <C is concerned. 

To get at this low-energy physics^jt, is convenient to 
use a strong-disorder RG approacbBti'Ll in which we it- 
eratively diagonalize the most strongly coupled parts of 
the Hamiltonian: Thus, at each step, we focus on the 
strongest bond |</ max | = ^ hi the system. Consider, 
first, the case when this is an antiferromagnetic bond: 
If this bond connects two spin-l/2s (such a bond can 
clearly be either even or odd), we freeze the two spin- 
l/2s into a singlet state, and introduce a renormalized 
coupling between the neighboring spins on either side: 
J = JlJrI Jmax.] here Jl and Jr are the bonds immedi- 
ately adjacent to J max , respectively to its left and to its 
right. If this bond connects two spin-Is (in this case, </ max 
must necessarily be an even bond), we again freeze them 
into their singlet ground state, and couple the neighbor- 
ing spins on either side to each other with the renormal- 
ized bond J = JlJr/ Jmax- Finally, consider the case 
where J m ax couples a spin-1 with a spin-1/2 (again, in 
this case, </ max must be even), with this spin-1/2 coupled 
to its other neighbor with a bond Jhaif, and this spin-1 
coupled to its other neighbor with a bond J on e- In this 
case, we form a renormalized spin-1/2 object representing 
the doublet ground state of this pair. The neighboring 
bonds Jhaif and J one remain unchanged in magnitude af- 
ter the RG step, but Jhaif changes sign. Finally, if J max 
is ferromagnetic (in this case, J max has to be an odd 
bond connecting two spin-l/2s), we put this pair of spin- 
l/2s into their triplet ground state, and replace it by an 
equivalent spin-1. 

This formulation of the RG is completely equivalent to 
the one used in Ref ||, and ignores O(l) coefficients in the 
RG recursion relations as well as the distinction between 
the bond J coupling a pair, and the gap in its spectrum — 
as in that case, I expect it to give accurate results when 
the effective value of disorder in the low-energy theory is 
large. My calculations will therefore be essentially exact 
in the regime of strong Griffiths effects (with a large value 
for at least one of the dynamical exponents), and are 
expected to be accurate even for smaller values of the 
dynamical exponents.Ej In the next two sections, I use 
this RG approach to analyze the low-energy physics of 
the effective model. 
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IV. RG FLOW-EQUATIONS AND FIXED-POINTS Finally, the last equation governs the T dependence of 



It is straightforward, if a little tedious, to write down 
the flow equations that describe the iterated action of the 
RG rules on our effective model, and this is what I turn 
to next. 

But first, some additional notation: Let us introduce 
Nff(^i), Nif([i), Nn(fi), and N a to denote the number 
of mass-/i //, If, II clusters, and a-type clusters respec- 
tively. Thus, we have 

N ff {n) = (1 - /„)(1 - g) 2 (l - q)q»N w , 
N lf Qt) = 2(1 - f a )g(l - g)(l - q)q»N w , 
N n (fl) = (l-f a )g 2 (l-q)q»N w , 

N a = f a N w . (3) 

In addition, let us denote the total number of //, If, 
and 11 clusters by Nff, N±f, and AT n respectively. It is 
also convenient to introduce the log-cutoff T = ln(f2o/f2), 
and the log-couplings (3 = ln(Sl/| J|), and keep track of 
the probability distributions of bond-strengths in terms 
of the probability distributions P s /o/w{P\^) of the cor- 
responding log-couplings. And finally, let us introduce 
P s/o/w( T ) to refer t0 P s/™/ o (0|r) respectively. 
The flow equations for the probabilities read 

^ =(1 - / Q )[P°(l - q)(l 5 2 (1 - /„)) - / a (P° + P°)] , 

^ =(1 - g)[P°(l - / a )(l - ,g 2 (l - g)) - g(P o + P s )] , 



while those for the functions P s / w / read: 

dr ~ dp + s s + 

+(qPo + (i - - fa)g 2 P°)(P s ® p s - P.), 

BP BP 

Ur ™ _ ur ™ , pOp , 

BY ~ 3f3 + ^ w + 

+(f a P° + (1 - g)(l - fa)g 2 P° s ){P w ®P W - P w ), 



(4) 



dP Q BP 



PoPo + 



BT 8(3 
+ (AP S ° + BP w ){P o ® P - P D ), 

where P ® P stands for the convolution 



(5) 



P®P = 



dp 1 d/3 2 P(/3 1 \T)P(f3 2 \T)S(p -Pi -02 



and A(T) and P(r) are defined as 

(l-/ Q )(l-.g(l-g)) 2 



A(T) = 
B(T) = 



2(l-g) + (2g-l)(f a + q)-2gqf a ' 

(l-g)(l-ff(l-/ Q )) 2 
2(l-g) + (2g-l)(f a + q)-2gqf a 



dN w 

dr 



-[i* + (1 - g)(l - /a)ff 2 P s U + /aP U ]^. (6) 



Some comments are in order before we go any further 
with our analysis: First of all, the very fact that these 
equations can be written down at all relies on the fact 
that of our ansatz, Eqn (JJj) for the mass distribution and 
relative abundances of various clusters, is consistent with 
the action of the RG transformation. Another thing to 
notice is that Eqns are invariant under the simulta- 
neous interchange of P, and P w on the one hand, and f a 
and q on the other. This is to be expected, and is related 
to the S — > —5 duality of the original problem. Indeed, if 
we consider interchanging P s and P w , which would corre- 
spond to the s type bonds being chosen weaker than the 
w type bonds, then our definition of the clusters would no 
longer be appropriate. In this case, it is more appropriate 
to think in terms of the complementary clusters enclosing 
w-type even bonds, and linked together by s-type even 
bonds. In this complementary description, the roles of 
f a and q will clearly be interchanged. Nevertheless, the 
(5^—<5 duality implies that the equations governing the 
flows in these complementary variables should have the 
same form as our earlier equations — in other words, the 
original equations must be form-invariant under a simul- 
taneous interchange of P s and P w , and f a and q. In the 
rest of the article, I therefore conform without any loss 
of generality to this suggestive labeling, and restrict at- 
tention to situations in which the s-type even bonds are 
on average stronger than the ui-type even bonds, and S 
is positive. 

The next order of business is clearly to check our ef- 
fective model and equations against the results of Ref ^ 
for the special S = case. Let us begin by setting 
P s = P w = P e , since there is no distinction at all between 
the u;-type and s-type even bonds at 5 — 0. Next, note 
that the average number of s-type even bonds is given 

by (1 - /„)(1 - q)N w g ( M + 1)5" = N w (l - / Q )/(l - q). 

Requiring that this equal N w gives the expected result: 
f a = q at S — 0. This immediately implies that the 
total number of antiferromagnetically coupled odd pairs 
of spin-l/2s equals 2f a N w . Furthermore, the number of 
ferromagnetically coupled odd pairs is expected to equal 
that of the antiferromagnetically coupled odd pairs on 
average. This implies that f a = {2Nff + N\f)/2N W — 
(1 — / a )(l — .<?)• Moreover, the ratio N of the number of 
spin- Is to the number of even bonds in the system (the 
notation is chosen to conform with that of Ref ^|) can 
be expressed as: N = (2N n + N lf )/2N W = (1 - f a )g. 
Putting all this together, we can express f a , g and q in 
terms of the single parameter N: f a = q=(l — N)/2, 
g = 2N/(1 + N). Furthermore, it is easy to check that 
Eqns (||,||) at 8 = are consistent with these relations, 
and reduce to three independent equations, two for the 
two functions P e and P , and one for the ratio N. In 
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fact, a little algebra shows that these reduced equations 
are completely equivalent to the equations of Ref [| 

The fixed points found in Ref ^ can now be translated 
into our language as follows: To begin with, the RSi 
phase boundary at R > R C (S = 0) is controlled by: 



9 = fa = (1 - AO/2 
g = 2N/Q. + N) = 1, 



0, 



P (fi\T) = Qe~^ 



P W ((3\T) = -e 



-/3/r 



(7) 



where Q is an arbitrary constant. [Of course, this is not 
really a fixed point for the probability distributions P s 
and P w , but becomes one if we transform to a descrip- 
tion in terms of the distributions of the scaled variable 
£ = /3/r; I will be sloppy about such terminology here 
and below, since it is always clear what is meant fr<j 
the context]. At such a fixed point, Eqn @ correctlj 
predicts N w ~ T~ 2 . 

On the other hand, the fixed points describing the 5 
Gapless Haldane phase are written as: 

q = f a = (1 - AO/2 = 1/2, 
g = 2N/(1 + N) = 0, 
p s (/?|r) = P W (P\T) - Per p (\ 
P o (0\T) = Qe-^ ,Q = 0, 



(8) 



with P an arbitrary constant; at any such fixed point, 
Eqn (^) implies My, ~ e _r , which is consistent with the 
previous results. EfQ 

Finally, the multicritical point at R = R c ,5 = is 
controlled by 

q = /„ = (1 - AT)/2 = 1/4, 
.g = 2N/(1 + N) = 2/3, 



P>(/?|r) 



P w (/3|r) = L -e^l\ 



(9) 



with Eqn (|6|) again correctlyBEl predicting A^ ~ r~ 3 . 

Let us now look for new fixed points that would cor- 
respond to the rather unusual Griffiths phases we expect 
to find in the general 5 > case. Consider first the 
(1,1) Griffiths phase at 8 ^ 0: In this phase, we ex- 
pect the odd bonds in the effective model to be much 
weaker than both types of even bonds (this would guar- 
antee that the resulting pattern of singlet bonds at long 
length-scales would have the topological order character- 
istic of the (1,1) phase). Moreover, we expect the typical 
mass of clusters to scale to zero at low-energies. Guided 
by such considerations, it is easy to see that our equa- 
tions admit the following two-parameter family of fixed 
point solutions that have all the 'right' properties: 

fa = 1, 



9 = 5 = 0, 

P s (/3|r) = P s e~ p »P, 
P W (P\T) = P w e~ p ^, 
Po(P\T) = Qe- Ql3 , Q 



= 0, 



(10) 



where P s and P w are two otherwise arbitrary positive 
constants that obey P s > P w (this condition matters for 
the stability of these fixed points). [Note that using these 
fixed point values in Eqn (^) implies N w ~ e~ p, " r .] 

Similarly, in the (2,0) Griffiths phase, we expect the 
u>-type even bonds to be much weaker than either the s- 
type even bonds or the odd bonds. In addition, we again 
expect the masses of clusters to go to zero at low-energies. 
Furthermore, from the physical picture developed earlier, 
we expect that there are two distinct regions in this Grif- 
fiths phase, a spin-1/2 rich region closer to the RS1/2 
phase boundary, and a spin-1 rich region adjacent to the 
RSi phase boundary. Again, it is not hard to find the cor- 
responding two-parameter families of fixed points. The 
region closer to the RS]^ phase boundary is controlled 
by: 



fa = 9 

q = 0, 

win 

Po(P\T) 
P„(0\D 



= 1, 



P s e- P ^, 
Poe-™, 
Qe- Ql3 , Q 



0, 



(11) 



where P s and P are otherwise arbitrary positive con- 
stants satisfying P s > P (again this controls the stabil- 
ity of these fixed points). [Note that using these fixed 
point values in Eqn (^) implies N w ~ e~ p ° r .] 

On the other hand, the region closer to the RSi phase 
boundary is controlled by: 



fa = q 
9 = 1, 

Po(P\T) 
PJ0\T) 



= 0, 



P s e- P °^ 
Poe-™, 
Qe~ QI3 , Q 



= 0, 



(12) 



with P s < P Q in this case (as before, this is related to 
the stability of these fixed points). [Note that using these 
fixed point values in Eqn (j6|) implies ~ e~ PsI \] 

Furthermore, for the 'degenerate' crossover case char- 
acterized by P s — P , we have an additional degree of 
freedom in the choice of f a . Of course, for any physical 
system, both f a and the common value, P, of P s and 
P w will be determined by complicated physics at higher 
energies — the corresponding functional relationship be- 
tween the two along the cross-over line in (R , 6) plane 
is clearly outside the scope of our analysis, although we 
will see later that it is possible to predict the leading 
behavior of both f a and P s as one approaches the mul- 
ticritical point along the crossover line. [Also, note that 
the non-universal value of f a along this line of fixed points 



7 



does not affect the scaling of N w , with Eqn (0) predicting 

Finally, the RS1/2 phase boundary is described by the 
last specimen in our menagerie of fixed points: 

fa = 1, 
q = 0, 

9 = 1/2, 

P s (f3\T) = P s e- p ^, 



P w (m = P„(J3\T) = -c 



-0/r 



(13) 



with P s an arbitrary positive constant. [Note that using 
these fixed point values in Eqn (|^) implies ~ r~ 2 , as 
expected.] 



V. SCALING FLOWS AND CROSSOVERS 

Naturally, the low-energy physics in any regime de- 
pends on the fixed point controlling it via the scaling 
flows and crossovers in the vicinity of the fixed point, 
and it is crucial to characterize these in order to develop 
a complete picture of the low-energy behavior. 



A. Near the Multicritical point 

Let us begin in the neighborhood of the multicritical 
point. From the (5 = analysis of Ref ^, we know that 
the flows away from this point along the R axis are con- 
trolled by the relevant eigenvalue A r — (\/l3 — l)/2. A 
more general stability analysis should yield another rele- 
vant eigenvalue A5, reflecting the expected instability of 
this fixed point to infinitesimal dimerization. Together, 
these two exponents will then control the multicritical 
scaling near this point in the usual way; in particular, 
the ratio A r /Aa will determine the shape of the RS1/2 
phase boundaries asymptotically close to the multicriti- 
cal point. 

Following Refs |||| let us look for eigenperturbations 
of the form 



p.(/?|r) = p s fp (/?|r) + 2( £e + estu )r^ 1 (i-^; 



-2/3/r 



P W (J3\T) = P:P(/3|r) + 2(6 e - £stu )r A - 1 (l-^)e-W r , 

r 

fa = fX + (e aq - e„/4)r\ 
q = q tp - (e aq + e„/4)r\ 

9 = g f P + (4e„/9 + 2e 3 /3)r\ (14) 

where the superscript fp denotes the values at the multi- 
critical fixed point. In the above, the coefficients e have 
been chosen so as to separate out the effects of going 



P a {fi\T) = pfp( / 3|r) + 2 £o r A - 1 (i--£ ) e 



away from 6 = from the effects of changing R at 5 = 0; 
thus, e aq = ( g = e sw = if 5 remains zero. Conversely, 
e„ and e e are both expected to be zero if R remains at its 
multicritical value. e is the only coefficient that can, in 
principle, get contributions both from deviations in R at 
(5 = 0, and from deviations in 6 for R fixed at its multi- 
critical value. However, one expect on physical grounds 
that the distributions of the odd bonds, related as it is 
to the effective couplings between the end-spins of a long 
segment in the Haldanc phase, is relatively insensitive to 
(5 so long as it is small — it is therefore reasonable to as- 
sume that e G is zero to linear order in S at fixed R. [In 
other words, the coefficients e are related to the pertur- 
bations in (5 and R as e = c r + . . ., e e = — c e r + . . ., and 
e S iu = c sw S + e aq = c aq 5 . . . and e„ = c n r + . . ., with 
all the c being positive and 0(1), and r = R— R c (0).] 

This ansatz is readily seen to solve the flow equations 
to linear order in the e, provided the following system of 
linear equations is satisfied: 



Ae„ 

Ae e 
Ae Q 

Ae ag 
Ae SU j 



-5e 9 , 

1 _ 3 _ 3 ^ 3 
6e g 2 6e 2 C °' 

1 1 1 

£g ~\~ £'n ) 

3 9 
2 f 

1 



= 2e a , 



(15) 



Fortunately, all eigenvalues of this system can be eas- 
ily determined as follows: To begin with, we can set 
e g = e aq = Csu, = 0, to obtain the three eigenvalues 
-1, -(1 + VT3)/2, and \ r = (VIE - l)/2 that charac- 
terize the (5 = flows; the corresponding eigenvectors 
live entirely in the (e e ,e ,e„), and correspond precisely 
to the eigenperturbations described in Ref ^. Next, note 
that A = —5 is clearly an eigenvalue; the correspond- 



ing eigenvector has e aq 



and lives entirely in 



the (e g , e n , e e , e Q ) subspace. Finally, we have two eigen- 
vectors that live entirely in the (e ag ,e SMJ ) subspace, with 
eigenvalues -(1 + \/l3)/2 and X 6 = (a/13 - l)/2. Thus, 
we indeed have just two relevant eigenvalues, the first A r 
corresponding to moving along the R axis at 6 = 0, and 
the other A5 representing the effects of non-zero dimeriza- 
tion at R = R c (5 = 0). The important, and at first sight 
rather surprising, feature is that these two eigenvalues 
are equal! In actual fact, this equality reflects an un- 
derlying 53 (permutation group of three elements) sym- 
metry of the multicritical point (corresponding to free 
interchange between the three phases that meet at this 
point); a discussion of this point, as well as generaliza- 
tions to some closely related problems, is the subject of 
a separate article.^ 

Now, a reasonable picture of the full crossover from the 
multicritical point to any of the phases in its vicinity is to 
assume that a system very close to the multicritical point 
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follows the multicritical flows as the RG proceeds until 
either rT Xr or5T x <> becomes 0(1), after which the system 
obeys the scaling flows characteristic of the fixed point 
governing the phase it is in. This matching procedure 
can then be used to relate the scaling behavior of various 
quantities near the multicritical point to the two relevant 
eigenvalues A r and Xs, and I will have more to say along 
these lines after we analyze the low-energy behavior of 
the different phases. For now, let us ask a more basic 
question, namely, the shape of the RS1/2 phase bound- 
aries close to the multicritical point. It is convenient to 
formulate the argument in terms of the scaling behavior 
of the some physical property, such as the topological or- 
der parameteM T = limr^oo T(r, 6, T) that characterizes 
the connectivity property of the (1,1) phase (see Ref 
for the precise definition). Multicritical scaling, in con- 
junction with the results of Ref [?] for 6 = 0, implies that 
we may write 



6~ 



T(r, 6,T) = T-^-^^FirT^, ) 



r 



(16) 



,F(x,0) 



6-24, 



x\ x -r , F(0, 0) a constant, 
m = -\/5. Now, if we keep r fixed at some small 



with lim 
and 

negative value and increase 8 from zero, the presence 
of the RS1/2 phase boundary at which the topologi- 
cal order is lost must be reflected in the the large T 
limit of F. Indeed, as we will see in the next section, 
one expects T(T) to scale as r-( 4 - 2<te ) on the RS 1/2 
phase boundary, with <f>jis = (y/E + l)/2. This can 
be consistent with the multicritical scaling form F only 
if lim^^-oo F(x, —y c ) ~ x (3-V5)/\ r f or sorne particular 

y c > 0, implying that <5 A <s = — y c r is the equation for the 
phase boundary close to the multicritical point. 
Since we have X r — Xs, this means that the phase bound- 
ary comes in linearly as shown in the schematic phase dia- 
gram Fig [l]. Of course, Fig |l| contains another ingredient, 
namely that the RSi/ 2 line slopes downwards, ruling out 
reentrant behavior. To justify this needs a somewhat 
more detailed consideration of the actual crossovers — 
while it is possible to do this using eigenvector infor- 
mation and the matching procedure outlined above, I do 
not pursue this further here. 



B. Near the RS1/2 and RSi phase boundaries 

The analysis in the vicinity of the RS^ and RSi fixed 
lines is much simpler, and the two cases are closely anal- 
ogous. Here, I discuss only the RS1/2 phase boundary in 
detail, confining myself to a brief summary of the corre- 
sponding results in the RSi case. 

In order to study small perturbations around a point 
on the RS1/2 line, it is convenient to parametrize these 
deviations as 

P s (/3|r) = P s e- p ^[l + 6.(1 - P s p)}, 



P w (P\r) 

Po(0\T) 



-PIT a 
[l + (S P -6 ow )(l-^)}, 



r 

e -P/T 



r 



[l + (6 P + 6 ow )(l-£)}, 



1 



9=^ + 5 g ), 

fa = l~ S a , 
q = 6 q . 



(17) 



Requiring that this ansatz satisfy the linearized flow 
equations yields a linear system of ordinary differential 
equations for the functions 6(T). Fortunately, the equa- 
tions for 5 a and 6 q are particularly simple: 



dT 
dS q 



4r ' 



(18) 



These immediately fix the T dependence of 8 a and 6 q to 
be 



6 a = C a T e 



-P S T 



S 9 = ( 



c a r 2 



r ' 



(19) 



at large T, both 5 a and 6 q thus decay extremely rapidly 
to zero (compared to the T x behavior we anticipate for 
8 w, S p , and 5 g ). Furthermore, it is easy to check that 
the linearized equation for S s only involves S q and 6 a : 



d5 s 
dT 



<Jq 

r 



6a_ 

4r' 



(20) 



which implies that the leading large T behavior of S s is 



6 S (T)~ const. + f^e- Psr . 



(21) 



The presence of the constant term in the above merely 
reflects the fact that we have a whole line of RS1/2 fixed 
points with P s (0|r) varying continuously along the line 
(in RG terms, this is a marginal coupling); requiring that 
the fixed point value of P s (0|T) be P 8 at the particular 
point about which we are perturbing tunes the constant 
to zero, implying that 8 S also decays very rapidly at large 

r. 

This simplifies matters considerably — for clearly, we 
are at liberty to set 6 a , 6 Sl and 5 q to zero in the linearized 
equations for 8 ow , 8 P , and 6 g in order to determine the 
eigenvalues A that control the slow growth or decay of 
these perturbations. The corresponding linearized equa- 
tions (with the 5 a / q / s set to zero) read: 



^ - -(26 

dr ~ r [ do 



2<U 



d5 p 
~dT 

~d~F 



ow 

r~ 



(22) 
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As in the previous section, the eigenmodes are clearly 
of the form 8 g / p / ow = e 5 / p / oto r A , with the coefficients e 
constrained by 

Xc ow — ~t~^OU'; 

Xe g = -2e g + 2e ow . (23) 

The three eigenvalues are +1, —1, and —2, with the 
respective eigenvectors proportional to [e p = , e ow — 
1, e g = 2/3), (e p = 1, e ow = 0, e g = 0), and (e p = 
0, e ow = 0, e g = 1). 

Thus, as expected, there is a single relevant eigenper- 
turbation on the RS^ line. Clearly, this corresponds 
to tuning the bare dimerization away from the critical 
value (5 = <5 C ; from the form of the relevant eigenvector, 
it is clear that this gives rise to a slight imbalance in the 
strengths of the w-type even bonds and the odd bonds 
in the effective model, which in turn drives g away from 
its critical value of 1/2 (as it must, since the fixed point 
value of g is 1 in the (2,0) phase, and in the (1,1) phase). 
This analysis also allows us to identify the appropriate 
crossover energy scale for a system close to the RSi/2 
phase boundary — clearly, any such system will 'look' 
critical for T less than the log-energy scale Ts eH ~ SZa, 
and will have behavior characteristic of either the (2,0) or 
(1,1) phase for T greater than this crossover scale (here 
S c ff =6 — S c )). All of the above is clearly-identical to the 
original analysis of the spin-1/2 chainJ§L3 and thus, the 
analysis here confirms the heuristic picture of Section [n|. 
Moreover, this also makes it clear that scaling behavior 
of the topological order parameter T in the vicinity of 
the RS1/2 phase boundary will be identical to the behav- 
ior of the dimerization order parametert-3 in the random 
spin-1/2 HAF chain in the vicinity of its RS1/2 critical 
point at zero dimerization; in particular, one therefore 
expect that T{R = R C ,S = S C ,T) ~ r( 4 " 2 ^ s )pn our 
RS1/2 phase boundary, with 4>rs = (\/5+ l)/2£3 

The analysis near the RSi line proceeds similarly: 
From the linearized flow equations, it is easy to establish 
that perturbations in the values of g, f a , and q all die 
away exponentially (apart from some unimportant pref- 
actors of powers of T) at large T. Furthermore, P ((3\T) 
also settles down to its fixed point form just as rapidly 
(of course, P o (0|r) is now a marginal coupling along this 
line, and the particular value P Q it settles down to de- 
pends on the initial condition, analogous to the behavior 
of P s in the RS]^ case above). This means we are at 
liberty to set all of these variables to their fixed point 
values in our analysis of the much slower (~ T x ) growth 
or decay of perturbations of P s and P w away from their 
common fixed point form. Using the usual parametriza- 
tion, 

PsW) = ~e-£(l + (e e + e sw )T\l £)), 
P W W) = ~e-£(l + (e e - e sw )T\l ~ £)), (24) 



for these slow deviations makes it clear that we again 
have precisely one relevant eigenvalue A = 1, with the 
corresponding eigenvector proportional to (e e = , e sw = 
1) (with the irrelevant eigenvector, of eigenvalue A = — 1, 
being proportional to (e e = 1, e sw = 0)). Moreover, 
the relevant perturbation clearly corresponds to turning 
on a slight dimerization, which introduces an imbalance 
in the strengths of the s-type and w-type even bonds. 
Finally, the corresponding crossover energy scale T$ again 
behaves has Tg ~ for small 6, in complete analogy 
with the spin-1/2 case. 

C. In the (1,1) Griffiths phase 

As we have seen earlier, the (1,1) Griffiths phase for 
S > is described by a two-parameter family of fixed 
points, with P s = P s (0\T and P w = P w (0\T) allowed to 
vary independently of each other, modulo the constraint 
-Ps > Pw As expected, a formal stability analysis at 
any point with P s > P w yields two marginal eigenper- 
turbations (with A = 0) corresponding to these two free 
parameters, with-all other perturbations dying away ex- 
ponentially in rJ13 However, most low energy properties 
in the phase are controlled by precisely these other terms 
that decay exponentially in T — after all, the exponen- 
tial decay in the log-energy T translates to power-law 
behavior in the energy f2, and it is precisely these power- 
laws that are the characteristic signature of any Griffiths 
phase. 

To analyze the low-energy flows in sufficient detail to 
get at this behavior, it is useful to parametrize small 
perturbations away from any S > fixed point labeled as 
{P S ,P W ) (with P s > P w ) as 

P S ((3\T) = P s e- p ^(l + e s (l - P.0)), 
P w (j3\T) = P w e- p -?(l + e w (l - P w 0)), 
P (f3\T) - P (T)e- p °^(l + e (l - P (T)f3)), 

9 = £g, 
fa 1 e ai 

Q = e q , (25) 

with -P (r) = C e~ p ™ r — this choice of P (T) 'factors' out 
the r dependence expected 'at' the fixed point, with e 
representing a small sub-dominant correction. Moreover, 
this ansatz clearly satisfies the flow equations to leading 
order in the e so long as the e obey a corresponding sys- 
tem of ordinary differential equations that governs their 
r dependence. 

Now, a simple analysis of the form of these equations 
reveals that the perturbations e s , e w , and e and their T 
dependence play no role in determining the leading large 
r asymptotics of e ff , e a , and e q . The e s / w / a can thus 
be set to zero as far as the analysis of the other pertur- 
bations is concerned, and the results of this analysis for 
e g/a/q can then be 'fed back in' to work out the leading 
behavior of the e s / w / (this latter step is actually quite 
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unimportant, for it turns out that the low-energy proper- 
ties of the system depend crucially on the behavior of the 
tg/a/q, and very little on that of e s / w / a ). This simplifies 
the equations for e g / a /q considerably, yielding: 



de a 
dT 
ckg_ 
dT 
de q 



(26) 



where, on the right hand side of each equation, only 
the terms that matter the most for the large T asymp- 
totics have been kept. The solution to these is relatively 
straightforward to write down: 



e g = (TC + C g )e- p - T , 



Ca{e -(p,-p w) r_ Cqe -P B r )j 



(27) 



where C Q , C q , and C g are constants of integration that 
depend on the initial conditions. Translating from log- 
energies to the energy scale f2 and keeping only the lead- 
ing power-law in each case, we thus obtain 



e s ~(C s + C>(fi /O))fi J 



n 



9 

P.-Pw 



(28) 



With this in hand, one can now quite easily work out the 
leading T dependence of e s / w / ; however, as this will play 
no role in our later discussion, I refrain from displaying 
the results of this analysis. 

To summarize, Eqn |2^ implies a rather simple picture 
for the effective Hamiltonian at cutoff scale ft: This pic- 
ture is in terms of type-a clusters, and 0-mass clusters 
of the 11, 1/, and // types, all connected to each other 
by w-type even bonds. Using the leading O dependence 
of N W) N w ~ Pro , derived earlier, the abundances of 
various types of clusters are seen to scale as 



N a 
Nff 



n 1 



•vP, +P„ 



N lf ~ (c H-£) + c g )n* 

JVn ~ «7> 2 (^) + 2C 9 C7 ( ,ln(^) + rt)<>'" (2!)) 



while the magnitudes of the exchange couplings all obey 
power-law distributions P(\ J\) ~ | J\~ 1+x , with x = P s 
for s-type even bonds, x = P w for w-type even bonds, and 
x ~ Q, Pw for the odd bonds. The picture that emerges is 
thus very similar to the results of the heuristic argument 
in Section [n|, with the exponents z X j 2 and Z\ introduced 
there given in terms of P s and P w as z\j 2 = , and 
z\ = (P s + Pw)^ 1 (in making this identification, I am 



of course ignoring the multiplicative logarithmic correc- 
tions, as well as the sub-dominant power-laws predicted 
by the more detailed analysis here). Furthermore, since 
Ps > Pw, Z\ is indeed less than z 1 / 2 /2, exactly as pre- 
dicted by the earlier Griffiths arguments. 

With this in hand, we can now match these low-energy 
flows with our earlier results for the critical and multi- 
critical flows to develop a reasonably accurate picture 
of the full crossovers as a function of f2 for a system 
in the (1,1) phase close to the multicritical point or the 
RSi/2 phase boundary. Consider first a system in the 
(1,1) phase, but very close to the RS^ phase boundary 
(and away from the multicritical point). In this case, the 
system will look critical for T less than the crossover value 
-TiScft ~ l^cffl 1 (with <5 ff = 5 — S c ), and then cross over to 
the flows described above at lower energies (i.e higher T). 
For instance, both P o (0|r) and Pn,(0|r) will scale down 
as r _1 for r above this crossover scale, while P s (0|r) 
will remain roughly constant. Beyond this point, P o (0|r) 
will decay rapidly, P o (0|r) - |<5 cff le" 1 ^ 1 ^" 1 ^ \ while 
P w / s (0\T) will both remain roughly fixed at their values 
at the crossover scale. Thus, the continuously varying 
exponent z 1 / 2 does indeed diverge as we approach the 
RS1/2 phase boundary, scaling as 5~g, while Z\ is smooth 
across the phase boundary. Moreover, applying a simi- 
lar argument to the function T(R,S,T) yields the scal- 
ing behavior of the topological order parameter T close 
to the RS1/2 phase boundary: T(R, 5, F) is expected to 
scale as Y~^ 2 ^ RS ^ for T > Tg D{l , and remains roughly 
constant as the energy is lowered further — as a result, 
T = limr^oo T(R, 5, V) will scale as TJ^~ 2lpRs) implying 



\5 e fi 



(4-20 RS ) 



(30) 



in the vicinity of the RS1/2 phase boundary. 

Next, consider a system in the (1,1) phase for S > 0, 
but very close to the multicritical point (the behav- 
ior at S = has already been discussed in Ref || and 
Ref 0). Matching the multicritical flows at higher ener- 
gies with the asymptotic behavior characteristic of the 
phase tells us that P s / w / o (0|r) all satisfy similar scal- 
ing forms: P s / w / o (0\T) = T- l K s / w / {rT~t , 5/r), where 
r = R — R c {5 = 0), and the second argument of the 
scaling function is determined by the fact that As = X r . 
Furthermore, requiring that this be consistent with the 
expected behavior in the (1,1) phase immediately implies 
that P s = lim r ^ooP(0|r) and P w = limr^oo P w (0\T) 
can be written as P s / W — \r\~Xr &~/ w ($ Vl r ll when r < 
(also note that the earlier results at S = 0,i'§ imply that 
$"^(0) must be a finite constant) 

In other words, the exponents Z\j 2 and z\ obey the 
scaling forms 



1 



Zl/2 = 



Zl 



ll/A. 



r| 1 0^1/2(^1) 



(31) 
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where Ej, 2 = l/$ w , H x = l/($ s + $„,), and r is as- 
sumed negative. Of course, in order to be consistent with 
our earlier analysis near the RS]^ phase boundary, the 
scaling function <I>~ must vanish as \y— y c \ when y = 6/\r\ 
approaches y c (corresponding to the RS]^ phase bound- 
ary), while $J must go smoothly to a constant. More- 
over, the present analysis tells us something else: Since 
this RS1/2 scaling sets in roughly when \y — y c \ < 1, 
the width in 6 of the critical regime controlled by the 
the RS1/2 fixed points vanishes linearly with r as one 
approaches the multicritical point. Furthermore, if we 
write Zi/2 ~ a r \ $ — &c \ 1 in the vicinity of the RS !/ 2 crit- 
ical line, then the scaling form derived above implies that 
the amplitude a r obeys a r ~ |r| 1_x ^ as we approach the 
multicritical point along the RS1/2 phase boundary. 

D. In the (2,0) Griffiths phase 

The situation in the (2,0) Griffiths phase is quite simi- 
lar. Except on the crossover line, a formal stability anal- 
ysis again yields two marginal eigenperturbations, corre- 
sponding to the freedom of choice of the two 'coordinates' 
P s and P that parametrize the corresponding family of 
fixed pointsEJ Of course, as in the (1,1) case, we need 
to go beyond such a formal linear stability analysis, and 
describe the flows in greater detail, to get at the low- 
energy physics. As before, our task is rendered easier by 
the fact that the deviations of the probability distribu- 
tions P s , P w , and P w from their fixed point forms play 
no role in such an analysis as far as the leading large 
r behavior of the other perturbations is concerned, nor 
do they matter much in our later calculations of various 
physical quantities; they may therefore be set to zero in 
our calculations. 

Below, I summarize the results of such an analysis sep- 
arately for each regime of the (2,0) phase. 

1. The spin- 1/2 rich regime 

It is convenient to parametrize the deviations of / Q , g, 
and q from their fixed point values as 

fa 1 f (ti 

.9=1- e s , 

q = e q . (32) 

Using this parametrization in the flow equations for f a , 
g, and q (in conjunction with the fixed point form for the 
distributions P s / w / ), it is easy to see that the e obey 

"^jT = — (^"s ~ Po)ta, 

= ~Po e g + Pgtq + C w e P °^ , 

§ = -{P„ + P Q )e q + 2e g e a C w e~ p ° r , (33) 



where only terms that play a role in determining the 
leading large T asymptotics of the e have been kept. [Note 
that I have used the fixed point dependence Pu,(0|r) = 
C w e~ p ° T of P TO (0|r) that follows immediately from the 
flow equation for P w upon using the fixed point values 
for all other parameters.] These equations immediately 
imply 

e q = C a C w (C w ln 2 (^) + 2C g ln(^) + C q ) e -^ +p ^ r , 

e a = C a e-^- p ^, (34) 

where C a , C q , and C g are all constants of integration. 

This gives a rather simple picture of the low-energy 
effective Hamiltonian at scale f2: The description is again 
entirely in terms of type-a clusters, and 0-mass clusters 
of the 11, If, and // types, all connected to each other 
by w-type even bonds. Using the leading fl dependence 
of N w , N w ~ f2 p °, derived earlier, the abundances of 
various types of clusters are readily seen to scale as 

N a ~n p °, 

N lf ~(c w M^-) + c g )n p ° +p ° 1 

Nff~ (^,ln 2 (^) + 2C 9 a,ln(^) + C^+^, (35) 

while the magnitudes of the exchange couplings all obey 
power-law distributions P(\ J\) ~ | J\~ 1+x , with x = P s 
for s-type even bonds, x ~ Q p ° for w-type even bonds, 
and x = P for the odd bonds. This is clearly very rem- 
iniscent of the heuristic picture of Section |fi|, with the 
exponents z 1 / 2 and z\ introduced there given in terms 
of P s and P as Z1/2 = P^/ 1 , and z\ = P^T 1 (in making 
this identification, I am of course ignoring the multiplica- 
tive logarithmic corrections, as well as the sub-dominant 
power-laws predicted by the more detailed analysis here). 

Thus, as expected, this regime of the (2,0) phase is dual 
at low-energies to the (1,1) phase that is separated from 
it by the RS]y2 phase boundary; we can pass from one to 
the other by interchanging the roles of the w-type even 
bonds and the odd bonds. Of course, in complete anal- 
ogy with our earlier analysis of the (1,1) phase, the ex- 
ponents z\ and Z\i<z satisfy the multicritical scaling form 
Eqn |3l], and again, Zi/% diverges as a r \6 — <5 C | _1 close to 
the RS1/2 phase boundary, with the amplitude scaling as 

a r ~ |r| 1- ~ for r = R — R c small enough. 

2. The spin-1 rich regime 

It is convenient to parametrize the deviations of f a , g, 
and q from their fixed point values as 



12 



fa 



(36) 



Using the fixed point form for the distributions P s / w / 
in the flow equations, it is easy to see that the e must 
satisfy 



dT ~ 9 ' " " 



-P.V 



de g 

dr 

de a 
dT 



-(P, + P s )e q + +2e g C w e 



-p s r 



-(P -P s )e a + 2P s e 



(37) 



where only terms that play a role in determining the 
leading large T asymptotics of the e have been kept. [Note 
that I have used the fixed point dependence P^Ojr) = 
C w e~ PsF of P w (0\T) that follows immediately from the 
flow equation for P w upon using the fixed point values 
for all other parameters.] These equations immediately 
imply 

e q = (ClC a ln 2 (^) + 2C w C g ln(^) + , 

^0 ^ , \ _p r 



{C w C a \n{^-) + C g )e 



9 v-»-»— v fi 

f - c e - {p °- Ps)T 



(38) 



where C a , C q , and C g are all constants of integration. 

Once again, this gives a rather simple picture of the 
low-energy effective Hamiltonian at scale ft: As before, 
the description is entirely in terms of type-a clusters, and 
0-mass clusters of the 11, 1/, and // types, all connected 
to each other by w-type even bonds. Using the leading 
f2 dependence of N w , N w ~ il Ps , derived earlier, the 
abundances of various types of clusters are readily seen 
to scale as 

N a ~n p °, 

Nff - (0 2 (^) + 2C g C wa H^) + C* g )Sl p >+ 2P °, 

N lf ~(c wa H^-) + c g )n p ° +p °, 



(39) 



where we have defined C wa = C w C a . The magnitudes of 
the exchange couplings all obey power-law distributions 
P(|J|) ~ |J|~ 1+a: , with x = P s for s-type even bonds, 
x ~ il p ° for w-type even bonds, and x — P for the 
odd bonds. And finally, as in the spin-1/2 rich regime, 
the exponents 21/2 and z\ introduced in Section || can 
be expressed in terms of P s and P a as z 1 / 2 = Po 1 , and 
z\ = P.T 1 (again ignoring multiplicative logarithmic cor- 
rections, as well as sub-dominant power-laws). 

We may now match our earlier results in the vicinity of 
the RSi phase boundary with this low-energy picture to 
obtain the scaling of z\ and Z\n close to the the RSi line. 



In a system close to the RSi line (but not close to the 
multicritical point), and for T < T$ ~ S^ 1 , P s / w (0\T) will 
both scale as T -1 , while P o (0|r) will stay constant. Be- 
yond this crossover scale, P s / o (0|r) will both stay roughly 
constant, while P^Olr) will fall off exponentially with 
increasing T. This immediately implies that Z1/2 will go 
smoothly to a constant as one approaches the RSi phase 
boundary, while Z\ will diverge as |<5| -1 — again, this is 
consistent with the Griffiths arguments in Section ||. 

The situation is somewhat different close to the mul- 
ticritical point, with 6^0: In this case, P s /w o (0|r) all 
obey the scaling form introduced earlier, P s / w /o{^W) = 
T~ 1 K s / w / p LrT'^r ,5/r). The requirement that the x — > 
+00 limitcil of the K s / w / (x ) y) be consistent with the 
behavior expected in this regime of the (2,0) phase im- 
mediately implies that P s = limr^co P s (0|r) and P = 
limr^oo P o (0|r) can be written as P s/o = r~$+ /w (5/r) 
when r > 0. In other words, the exponents z 1 / 2 and Z\ 
obey the scaling forms 



1 



21/2 



zi 



, r |l/A r "1/2 r - 



ll/A, 



(40) 



where H+, 2 = 1/$^, = l/$f, and r is assumed pos- 
itive. Of course, for this to be consistent with our ear- 
lier analysis near the RSi phase boundary, we must have 
$o(0) finite and non-zero, and $J (x) vanishing linearly 
with x for small x. This also tells us something more 
about the RSi scaling for small r: Clearly, the width 
in 5 of the RSi critical region vanishes linearly with r 
for small r; furthermore, if we write z\ ~ a r |£| near 

the RSi phase boundary, then the critical amplitude a r 
1 1 

scales as r A >- for small r. 



3. Along the crossover line 

The analysis of the crossover case is very similar to 
that in the two regimes on either side, with only some 
minor differences. As before, the deviations of the three 
distributions P s / w / from their fixed point form plays no 
role in determining the the decay of the perturbations of 
g and q from their fixed point values. Furthermore, it 
is easy to see that the leading large T behavior of these 
perturbations is also independent of precisely how /„ ap- 
proaches its non-universal fixed point value /. 

Parametrizing them as 



9 = 1 - e s 



<-</■ 



(41) 



it is easy to see that they obey the following system of 
equations 
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dEg_ 

dT 
dc 



-Pe n + Pe n 



fC w e 



-pr 



= -2Pe q + 2(1 - f)C w e 



-pr. 



(42) 



where P denotes the common fixed point value of 
P s/o (0|r), C w e~ pr is the fixed point value of P^Oir), 
and only kept the terms that matter the most at large T 
have been kept. The leading large T behavior of e g and 
e q immediately follows: 



e q = (f(l-f)C^ + 2C w C g (l 
(9 = {C a + fC w T)e- pr , 



f)T- 



C q )e 



-2Pr 



(43) 



where C g and C q are constants of integration. This can 
now be used to determine the manner in which P s / w / 
settle down to their fixed point values, as well as follow 
the approach of f a to its fixed point value. However, since 
none of this matters for the leading low-energy behavior 
of the densities of various clusters, I do not pursue this 
any further here. 

These results for e g and e q give a simple picture of 
the low-energy effective Hamiltonian at scale f2: As in 
all the other cases, the description is entirely in terms of 
type-a clusters, and 0-mass clusters of the 11, If, and 
// types, all connected to each other by w-type even 
bonds. Using the leading O dependence of N w , N w ~ 
£l p , derived earlier, the abundances of various types of 
clusters are readily seen to scale as 





-m p , 


N n - 




Nff r 


-(1-J 


N lf r 


^2(1- 



0(A 2 ln 2 (^) 



2AC g ln( 



«0x 



c 2 g )n 3P : 



c g )n 



2P 



(44) 



where the constant A is defined as A = fC w . [Of course, 
as one moves along the crossover line and approaches the 
multicritical point, / will tend to the universal multicrit- 
ical value of 1/4, and P will vanish with the exponent 

l/V] 



VI. PHYSICAL PROPERTIES 

This detailed picture of the effective Hamiltonian in 
various parts of the phase diagram can be used, in 
principle, to obtain corresponding information about 
the low-temperature (T) thermodynamics, and the low- 
frequency, low-temperature dynamics. Below, I focus on 
the low temperature specific heat and susceptibility, as 
these provide us with a simple, experimentally verifiable 
signature of the unusual crossover from spin- 1/2 to spin-1 
behavior as we tune across the crossover line in the (2,0) 
phase. [Only the Griffiths phases are discussed, since the 
behavior at the multicritical point and-.the RSi phase 
boundary has been analyzed in earlieiBQ work, and the 



properties of the RS1/2 phase boundary are identical to 
that of a spin- 1/2 chain in its Random Singlet state, an- 
alyzed earlier by Fisher B\ 

As in Ref [| the broad power-law distributions of the 
couplings in the effective Hamiltonian make it possible 
for us to treat the low-temperature thermodynamics in a 
relatively simple way. The basic idea is to run the RG till 
the cutoff n is reduced to be Qt ~ T, and recognize that 
the temperature T dominates over almost all exchange 
couplings in this renormalized problem with cutoff f2y. 
As a result, the leading low-temperature behavior of the 
system can be understood in terms of the thermodynam- 
ics of free spin-1 and spin-1/2 objects, whose densities 
are obtained from our earlier results for the abundances 
of various clusters in the effective Hamiltonian at scale 

n T . 

For instance, in the (1,1) Griffiths phase, the leading 
contribution to the low-temperature entropy will come 
from the entropy of the almost free spin-1/2 degrees of 
freedom that make up the type-a clusters. 



C'v = Tln(2) x 2 



dN a (Q T ) 
dT 



(45) 



where we have set ks to one. This implies that the lead- 
ing temperature dependence of the specific heat will be 

1 

a non-universal power-law of T: CV ~ T* 1 / 2 . These 
spin-1/2 degrees of freedom also provide the dominant 
contribution to the zero-field susceptibility: 



X 



AT 



x 2iV a (r! T ), 



(46) 



where g is the gyromagnetic ratio, and we have also set 
/is, the Bohr magneton, to 1. The leading temperature 

dependence of \ is thus x ~ 1/^ * 1/2 

Forming the 'Wilson-ratio' W = Tx/g 2 Cv, it is clear 
that the low-temperature limit of W is 



W = C1/2Z1/2 



(47) 



with the universal constant C\/2 — 1/(4 ln(2)) character- 
istic of a spin-1/2 system. Since the spin-1 degrees of 
freedom are always sub-dominant in this phase, this re- 
sult holds throughout, and the situation for general 5 ^ 
is thus not very diffecftnt from that predicted for the 6 = 
case in earlier workJSu 

The (2,0) Griffiths phase, on the other hand, presents 
some unusual possibilities: Clearly, the above analysis 
carries over unchanged to the spin-1/2 rich regime of the 
(2,0) phase, and the Wilson ratio takes the same form as 
in the (1,1) phase. However, in the spin-1 rich regime, 
the thermodynamics in the low-temperature limit is dom- 
inated by type-11 clusters, and therefore controlled by 
the exponent z\ that determines their abundance: 



C v = Tln(3) x 2 



dNnjOr) 
dT 



(48) 
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which implies that Cy ~ T*i in the low temperature 
limit. Similarly, the susceptibility x may be written as 



X=% x2JV 11 (fi r ), 



(49) 



which implies that x ~ V^ 1 Z1 m ^ ne l° w temperature 
limit. In other words, we again have power-law depen- 
dences similar to the spin-1/2 rich regime, with the Wil- 
son ratio 



W =C lZl 



(50) 



again proportional to the dominant dynamical expo- 
nent. However, the constant of proportionality is now 
completely different, and takes on the universal value 
2/(3 ln(3)) characteristic of a spin-1 system. Finally, 
along the crossover line, one again predicts similar power- 
law dependences controlled by the common value, z, of 
Z\ and Zxji-, but the Wilson ratio will now be completely 
non-universal, with W = cz, and c varying continuously 
along this line. Thus, the low temperature measurements 
of the Wilson ratio at different points in the (2,0) phase 
provide one way of getting at the unusual physics of this 
phase in which the system looks, at low temperatures, 
like a spin-1 one chain in one regime, and a spin-1/2 
chain in the other. However, since the low temperature 
thermodynamics is always controlled by the dominant 
dynamical exponent, a direct experimental handle on the 
second dynamical exponent at any particular point in the 
phase is lacking. Of course, the first corrections to the 
leading power-law behavior are also straightforward to 
calculate, and the experimental data could therefore be 
fit to such a more detailed formula to obtain the sub- 
dominant dynamical exponent at any particular point; 
unfortunately, this would probably not be a particularly 
compelling test of the theoretical picture presented here. 

I conclude with some speculations regarding experi- 
ments that could be directly sensitive to physics con- 
trolled by the second dynamical exponent: One might 
imagine that both exponents will leave an imprint on 
some dynamical property like the inelastic neutron scat- 
tering cross-section. However, reasoning as in Ref |22|, it 
is clear that the low-frequency intensity of the main fea- 
ture in the spectrum (at wavevector it) will again be con- 
trolled by the dominant exponent, with the contribution 
of the sub-dominant degrees of freedom only providing a 
small correction to this dominant low-frequency behav- 
ior. While it is possible that the low-frequency intensity 
in some other parts of the Brillioun zone might have some 
signature of these sub-dominant contributions, this is not 
at all obvious, and remains only a tantalizing possibil- 
ity for now. Another possibility relates to the fact that 
strong ferromagnetic correlations between some pairs of 
widely separated spins are present in the ground state 
in either Griffiths phase. Since the occurrence of such 
correlations is controlled primarily by the second, sub- 
dominant exponent, any signature of these in the static 



structure factor at low-temperatures could also give us 
a way of directly measuring this exponent. Finally, in a 
regime in which the spin-l/2s dominate, the most direct 
way of seeing the spin- Is would be to use some dynami- 
cal probe that is preferentially sensitive to spin-1 degrees 
of freedom, via a 'selection-rule' requiring that Am z = 2 
for the transitions induced (where m z is the z projection 
of the spin quantum number of a state), but again, this 
has not yet been backed up by any specific calculations. 



VII. OUTLOOK AND CONCLUSIONS 

Thus, the low-energy behavior of random antiferro- 
magnetic spin-1 chain systems presents a particularly 
interesting example of Griffiths effects if there is some 
dimerization in the exchange couplings on average. In- 
deed, as I have demonstrated above, Griffiths effects in 
such systems lead to a rather unusual Griffiths phase — 
in one region of such a phase, the system looks like a 
spin-1/2 chain at low energies, while it transforms itself 
into a spin-1 chain in another part of the same phase. 
Somewhat fortunately, this unusual low-energy behav- 
ior may be described in detail, and some concrete ex- 
perimentally verifiable signatures in the low tempera- 
ture behavior identified, via an 'almost-exact' analyti- 
cal strong-disorder renormalization group approach. The 
RG approach presented here yields a detailed picture for 
the effective Hamiltonian valid at low energies in vari- 
ous parts of the phase diagram, and an intriguing pos- 
sibility for future work is to use this picture to identify 
some compelling dynamical signatures of this interesting 
low-energy physics. Of course, Griffiths phases similar 
to those described here will exist in higher spin chains 
as well. Given the recent results of Ref 12, S = 3/2 



chains are of particular interest in this regard, especially 
since such Griffiths effects could be important even at 
relatively low values of randomness in these systems (the 
GriffithS|-ohases of S — 3/2 chains will be discussed in a 
separate£3 article) . Finally, it would be interesting to ask 
if similar disorder effects may exist in two-dimensional 
magnets with more than one topologically distinct phases 
for the pure system. 
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